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I. INTRODUCTION 


With the worldwide development of the radioactive ion 
beam facilities since the 1980s, the study of the exotic nu¬ 
clei becomes a hot frontier of nuclear physics HI. The exotic 
nuclei is a kind of nuclei with extreme neutron-proton ratio, 
which exhibit quite different features from the stable nuclei 
or the ones near stability line. As one of the representatives 
with the nuclear novel phenomena, the bubble nucleus is an 
exotic one characterized by the distinct central depressions 
of the matter distributions, namely the bubble-like structure. 
The study of bubble nuclei can be traced to the pioneering 
work of Wilson as early as in the 1940s ||2], in which the nu¬ 
cleus was assumed to be a thin spherical shell to explain the 
equally spaced nuclear single-particle levels. After that in the 
1970s, a variety of nuclear models, including the liquid drop 
model 13, the Thomas-Fermi model Q, and the Hartree-Fock 
method 0, have been applied to test the existence of the nu¬ 
clear bubble-like structure. Usually, it is recognized that the 
bubble-like structure originates from the low occupancy of s- 
orbit near the Fermi surface, which is the only wave function 
with non-zero value at r = 0. 

The nuclear charge distribution is an important observable 
which can provide very detailed information of nuclear struc¬ 
ture IMl- For instance, the charge distribution can reflect 
the proton-density distribution in nucleus, and a central de¬ 
pressed charge distribution is the consequence of the proton 
bubble-like structure. Hence, a proton-bubble nucleus can be 
identified experimentally from the measurement of the charge 
distribution, e.g., by the elastic electron-nucleus scattering ex¬ 
periment. In the near future, more experimental data about the 
charge distribution of exotic nuclei are expected to be obtained 
both from the project of the SCRIT ||9l and ELISe IfTOl . 

In medium mass region of the nuclear chart, some candi¬ 
dates of the proton-bubble nuclei were predicted, such as ^'^Si 
ifTTHTSl and "^®Ar ifTsUTTIl . For '^^Ar, the occurrence of proton 
bubble-like structure is found to essentially depend on the or¬ 
der of the proton orbits n 2 s\i 2 and while the prediction 

is evidently model dependent. In the calculations of the non- 
relativistic Skyrme Hartree-Fock-Bogoliubov(HFB) approach 
ini, it is found that the proton bubble-like structure may 
emerge in "^®Ar as the conclusion of the inversion of proton 
7r2si/2 and nld^ji orbits with tensor interaction included. Sim¬ 
ilarly, the relativistic mean-field theory (RMF) models with 
nonlinear meson self-couplings MM predict also a highly 
depressed density profile for ^®Ar without the tensor force. 
On the other hand, the pairing correlation effect could signif¬ 
icantly quench the bubble-like structure. In Ref. ca, with 
semirealistic NN interaction, the proton bubble-like structure 
is unlikely to exist in the argon isotopes due to the strong pair¬ 
ing effects. Different from ^®Ar, the emergence of the proton 
bubble-like structure in ^'^Si may be owing to the large gap be¬ 
tween the proton n2s\i2 and Jild^i 2 states, namely the Z = 14 
proton shell IfTSl . However, as predicted by the particle- 
number and angular-momentum projected generator coordi¬ 
nate method (GCM) based on the mean-field approaches, the 
dynamical correlation might strongly quench the bubble-like 
structures in both ^'^Si 031201 and "^®Ar 1211 . Therefore it is 


still an open question whether the proton bubble-like structure 
exists in "^®Ar or ^'^Si. Besides these two, the central charge 
density of is also predicted to be depressed distinctly in 
Ref 03, thus it is worthwhile to study the systematics of the 
charge density profiles along the N — 2% chains. 

As addressed above, the emergence of the proton bubble¬ 
like structure in the nuclei around “^^Ar is not only tightly 
related to the proton configurations near the Fermi surface, 
namely the position of the proton orbits n2s and nld and 
the gaps between, but also to the effects of the pairing 
and dynamical correlations. In recent years, as the natu¬ 
ral extension of the density-dependent relativistic Hartree- 
Fock (DDRHF) theory ||23|23, the relativistic Hartree-Fock- 
Bogoliubov (RHFB) theory with density-dependent meson- 
nucleon couplings was established to describe the weakly 
bound exotic nuclei 11241 . With the presence of the Fock terms 
in the mean-field channel, especially the inclusion of p-tensor 
couplings, substantial improvements have been achieved in 
the self-consistent descriptions of the nuclear shell struc¬ 
ture and the evolution of the ordinary and superheavy nuclei 
Il23ll25ll26l . the nuclear novel phenomena including the halo 
structures in cerium and carbon isotopes llZ7ll28ll . and the nu¬ 
clear spin-isospin excitation modes lEaHii]. 

Inspired by the above mentioned advantages, in this work 
we take the argon isotopes and N - 2% isotones as the 
candidates to explore the occurrence of the proton bubble¬ 
like structure within the RHFB theory, as compared to the 
relativistic Hartree-Bogoliubov (RHB) calculations ll32H33 . 
Based on the calculations with both RHFB and RHB theories, 
the bulk properties, the charge-density profiles, and the proton 
single-particle levels of the selected nuclei will be analyzed in 
details. The contents are organized as follows. In Sec. 
we briefly recall the general formalism of the RHFB theory 
and the charge-density profile. In Sec. Ill are given the de¬ 
tailed discussions on the emergence of the proton bubble-like 
structures in the selected isotopes and isotones, including the 
model-reliability test, the charge-density profiles, the proton 
single-particle spectra, and the relevant pseudo-spin symme¬ 


try. Finally, a brief summary is drawn in Sec. IV 


II. THE METHOD 


Under the framework of the DDRHF theory and its exten¬ 
sion — the RHFB theory, the nucleons in finite nuclei are 
treated as point-like particles moving independently in the 
mean field provided by the others, i.e., the mean-field ap¬ 
proach. Consistent with the meson-exchange picture of nu¬ 
clear force, the model Lagrangian as the theoretical starting 
point is then constructed by including the degrees of free¬ 
dom associated with the nucleon (ij/), the isoscalar cr- and lo- 
mesons, the isovector p- and ;7r-mesons, and the photon (A) 
12^133 . At the level of mean-field approach, the isoscalar 
cr- and cu-meson fields dominate the nuclear attractive and 
repulsive interactions, respectively, and the isovector p- and 
TT-meson degrees of freedom are responsible for the isospin- 
related aspects of nuclear force, and the photon field stands 
for the electro-magnetic interactions between protons. In fact. 
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not only the isovector ones, the Fock diagrams of the isoscalar 
cr- and w-couplings also carry the isovector nature of nuclear 
force, showing substantial contributions to the symmetry en¬ 
ergy OH [381. 

Following the standard variational procedure, one can de¬ 
rive the field equations of nucleon, mesons and photon from 
the Lagrangian, i.e., the Dirac, Klein-Gordon and Proca equa¬ 
tions, respectively. Meanwhile, the continuity equation, lead¬ 
ing to the energy-momentum conservation, can be also ob¬ 
tained, from which is derived the Hamiltonian of the system. 
In this work, which deals with the nuclear ground states, the 
time-component of the four-momentum carried by the mesons 
are dropped, which amounts to neglecting the retardation ef¬ 
fects in the Fock terms ll^ . Substituting the field equations 
of mesons and photon, the Hamiltonian can be expressed as 
the functional of nucleon field. In terms of the creation and 
annihilation operators (cj,, Ca) defined by the stationary solu¬ 
tions of the Dirac equation, the Hamiltonian operator can be 
further quantized as 

^ clclcp'Ca' Yj Kna'P” ( 1 ) 

oyS aa'PP' 0 

where Tap represents the kinetic energy term, and the two- 
body potential energy terms correspond to the meson- 

(or photon-) nucleon couplings denoted by (p. 

Tap = J dri//a(r)(-ir ■ V + M)iffp{r), (2) 

Kpa'p' drdr'ipa(r)ipp(r')r^(r,r') 

xD^(r,r')i/rp,(r')ipa>(r). (3) 

In the above expressions, ipa stands for the Dirac spinor of 
nucleon, M is the nucleon mass, T^ir,r') corresponds to the 
interaction matrices of various meson-nucleon coupling chan¬ 
nels, i.e., the cr-scalar, w-vector, p-vector, p-tensor, p-vector- 
tensor, 7r-pseudo-vector, and photon-vector couplings, and 
D^{r,r') represents the meson (photon) propagator. 

In the limit of mean-field approach, the contributions from 
the negative energy states are generally neglected, namely the 
no-sea approximation. The nuclear energy functional E is 
then determined as the expectation of the quantized Hamil¬ 
tonian Q with respect to the Hartree-Fock ground state 

E = (OolHIOo), |Oo) = Y] cJlO), (4) 

i 

in which the index / denotes the states corresponding to pos¬ 
itive energy and |0) is the vacuum state. In contrast to the 
RHB theory, the RHFB theory includes both the Haitree and 
Fock mean fields which correspond to the direct and exchange 
terms in the expectation of the two-body potential with re¬ 
spect to the ground state |Oo), respectively. 

For spherical nuclei, the relativistic Hartree-Fock (RHF) 
equation can be derived from the variation of the energy func¬ 
tional E 0 with respect to the Dirac spinor i/f(r) as d: 

I dr'h{r,r')il/{r') - e\J/{r), (5) 


where e is the single-particle energy and the single-particle 
Hamiltonian h{r, r') consists of three terms: the kinetic energy 
part the local and nonlocal potentials, i.e., and h^, 
respectively. 


r') =(ar ■ p H- pM)6(r - r'), 
h%,r') =[SH/-)r5 + SoCr) + /5;s(r)]5(r - r'), 
. , (YG(r,r') Yf(r,r')] 

XG(r,r ) XF(r,r ) j 


(6a) 

(6b) 

(6c) 


In /i^, the local potentials Sq and include the Haitree 
mean fields and the rearrangement term, and the nonlocal ones 
Yg, Yf, Xg and Xf in correspond to the Fock mean fields 

El- 

For the open shell nuclei, one has to take the pairing corre¬ 
lations into account. In order to provide a reliable description, 
we incorporate the Bogoliubov scheme with the DDRHF the¬ 
ory to deal with the pairing effects, leading to the RHFB the¬ 
ory 12411 . Following the standard procedure of the Bogoliubov 
transformation l^l40l . the RHFB equation can be derived as: 


f 


dr' 


' h{r,r') 
-A(r,r') 


A(r,r')'l 

h(r,r'). 


>ffu(r'y 


'A + Eq O' 

'^uiry 

1 

o 



(7) 


where ipG and [j/y represent the quasi-particle spinors, E^ is 
the single quasi-particle energy, the chemical potential A is 
introduced to preserve the particle number on the average, and 
the pairing potential A(r, r') reads as: 

A„(r,r') = - ^ y Vl^p{r,r')Kp{r,r'), (8) 

A 


with the pairing tensor k, 

^„(r,r')-(Av„(r)>f/„(r'). (9) 

In the particle-particle {pp) channel, we utilize the finite-range 
Gogny force DIS ST] as the effective pairing force. Aim¬ 
ing at the nuclei around ^®Ar, the original Gogny force DIS 
can provide appropriate descriptions of the pairing effects as 
demonstrated in Ref. sa. Notice that the RHFB equation 
Q is a coupled integro-differential equation and is hard to 
be solved in coordinate space. In order to provide an appro¬ 
priate description of the asymptotic behaviors of density pro¬ 
file, we expand the quasi-particle spinors tj/G and ij/v on the 
Dirac Woods-Saxon (DWS) basis f43l . and the basis param¬ 
eters, namely the spherical box size R^ax and the numbers of 
positive and negative energy states (resp. Nf and No), are 
taken as = 28 fm, Nf - 44, No — 12. 

In this work, the charge density is determined from the 
proton-density profile by incorporating the corrections of the 
center-of-mass motion and finite nucleon size. The first cor¬ 
rection is done by using the proton density in the center- 
of-mass reference frame, i.e., pc.m. which is related to the 
Hartree-Fock (HF) proton density through, 

Phf(0 = f r'Ddr', (10) 

B^7r2 J 
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where = 2{Pl^)l(3tP') and Pc.m. is the center-of-mass 
momentum. The second correction is taken into account by 
doing the convolution of Pc.m. with a Gaussian representing 
the form factor. 


Pch(r) = 


1 



k sin (kr)piif(k) exp 



dk, 


( 11 ) 

where Phf(^) is the Fourier transform of the HF proton density 
and — 2/3(0.862^ - 0.336^NIZ) accounts for the finite nu¬ 
cleon size nil. Denoting/!^ = l/(a^-B^), the charge-density 
distribution pch is finally derived as, 


Pch(r) = 


— ^ r r'dr'piip(r')\i 


g-F(r-r')- _ g-i"(r+r' 


( 12 ) 


where Phf(^ 0 corresponds to the proton density determined 
by the self-consistent calculations with the RH(F)B theories. 


III. RESULTS AND DICUSSIONS 


In this paper, we focus on the occurrence of central depres¬ 
sions of the charge-density profiles of the nuclei around “^^Ar, 
i.e., the proton bubble-like structure. The calculations are per¬ 
formed with the RHFB and RHB theories using the optimal 
effective interactions on the market, namely the RHF ones 
PKAl 1231, PKOl 1221 and PK03 El, and the RMF ones 
PKDD El and DD-ME2 l47l . 

Taking argon isotopes as the representatives, we firstly test 
the model reliability in terms of the binding energies and the 
root-mean-sqaure (rms) charge radii, as referred to the exper¬ 
imental data 1481 l49l . Figures[TJa) and[2b) display the devia¬ 
tions of the calculated binding energy per particle 1 A and 

charge radii from the data, respectively. It is found that 
both RHFB and RHB models provide appropriate agreements 
with the experimental data on the binding energies, whereas 
in the results of PKOl, PK03, PKDD and DD-ME2 notable 
discrepancies appear on the neutron-deficient side. Eor the 
charge radius which contains both corrections of center-of- 
mass motion and finite nucleon size, the DDRHE functional 
PKAl presents precise agreement with the data P9ll , and other 
selected models present relatively smaller values of rch than 
the data. Evidently, the RHEB theory with PKAl provides the 
most reliable descriptions on the bulk properties of the argon 
isotopes, particularly the charge radii. 

As we mentioned in the introduction, "^®Ar with two pro¬ 
tons deficient from '^^Ca was predicted as the candidate of a 
proton-bubble nucleus, which is characterized by the distinct 
central depression of charge-density distribution, if the inver¬ 
sion of the order of proton (tt) orbits n 2 s \/2 and nld-iji oc¬ 
curred iUMni. However, such inversion is essentially model 
dependent. To clarify the situation, we plot the charge-density 
profiles calculated by using the effective interactions PKAl, 
PKOl, PK03 and PKDD respectively in Pig. |^(a)-(d) for 
the argon isotopes. The results calculated with DD-ME2 is 
omitted, which shows similar systematics as PK03 [see Pig. 
|^(c)]. In Pig. 1^ it is clearly shown that the effective interac¬ 
tions PKOl, PK03 and PKDD, which seem to predict '^^Ar to 



14 16 18 20 22 24 26 28 

N 


FIG. 1. (Color online) Deviations of the calculated binding energies 
per particle £/A (MeV) [plot (a)] and charge radii tch (fm) [plot (b)] 
from the experimental data 148 1 1491 for argon isotopes. The theoret¬ 
ical results are extracted from the calculations of RHFB with PKAl, 
PKOl and PK03, and of RHB with PKDD and DD-ME2. 



r (fm) r (fm) 

FIG. 2. (Color online) The charge-density distributions of argon iso¬ 
topes calculated by RHFB with PKAl, PKOl, PK03 and RHB with 
PKDD. 
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have the bubble-like structure, present similar charge-density 
profiles with distinct central depressions for the nuclei around 
"^®Ar. On the contrary, the central depressions never appear 
in the charge-density profiles determined by PKAl along the 
argon isotopic chain. Notice the fact that PKAl provides the 
best agreement with the data of charge radii of argon isotopes 
among the selected models as shown in Fig. [3b). It seems 
that the occurrence of proton bubble-like structure is not fa¬ 
vored in the argon isotopes. 



FIG. 3. (Color online) Canonical proton single-particle energies for 
argon isotopes by RFIFB with PKAl and RHB with PKDD. The 
lengths of thick bars correspond to the occupation probabilities of 
proton orbits. 

To understand the charge-density profiles of the argon iso¬ 
topes, Fig. shows the proton canonical single-particle spec¬ 
tra calculated by RHFB with PKAl and RHB with PKDD, 
which present rather different charge distributions. In Fig. [3 
the lengths of the ultrathick bar denote the occupation proba¬ 
bilities of the orbits. As pointed out in Refs. lfT^[T3l[T7l . the 
order of the proton (tt) states 7rl£/3/2 and ;7r2ii/2 , as well as the 
gap between, is crucial for the occurrence of the bubble-like 
structure in "^®Ar. Along the isotopic chain of argon, it is found 
from Fig. [^a) that PKAl gives deeply bound and near fully 
occupied n 2 s\i 2 state, which does not support the formation of 
the bubble-like structure. While in Fig. j^b) the calculations 
with PKDD present an inversion on the order of the states 
7rl£/3/2 and Tils\j 2 at "^°Ar {N — 22) and after that the proton 
state 7r2ii/2 is less and less occupied, leading to the occur¬ 
rence of the central depressions of charge density in Ar [see 
Fig. |3d)]- Concerning the shell closures Z = 14 IITSl and 16 
lEol, PKAl presents distinct energy gap between 7rW3/2 and 


7r2si/2 that gives the shell Z - 16 at neutron-deficient side [see 
Fig. 13^)], and approaching the neutron-rich side, this shell 
(Z = 16) is strongly quenched and the one Z = 14 emerges, 
leading to well preserved pseudo-spin symmetry, i.e., nearly 
degenerated proton orbits n2s\i2 and nld 2 ,i 2 in "^®Ar (N = 28). 
On the contrary, the Z = 14 shell is persistent well in the cal¬ 
culations with PKDD along the isotopic chain. 

It is worthwhile to mention that the analysis of beta-decay 
Q values, single neutron separation energies, and the energies 
of the first excited state indicate the existence of the magic 
number Z - 16 in neutron-rich regions of nuclear chart ifSOl . 
To test the model reliability. Fig. [3 shows the proton single¬ 
particle energies of the sulfur isotopes from A = 16 to 30, 
calculated by RHFB with PKAl [Fig. [Ja)] and RHB with 
PKDD [Fig. |3!b)]. It is found that the calculations with PKAl 
give consistent prediction on the emergence of the proton shell 
Z = 16 with the analysis in Ref. ll50l . While similar as argon 
isotopes, the RMF Lagrangian PKDD predicts only the proton 
shell Z - 14 to occur in the proton spectra of sulfur isotopes. 
For the other selected Lagrangian, i.e., PKOl, PK03 and DD- 
ME2, similar proton spectra are predicted as PKDD for both 
argon and sulfur isotopes. 



FIG. 4. (Color online) Same as Fig. [3 but for sulfur isotopes. 

As discussed in Refs. nsHni, the occurrence of the cen¬ 
tral depression of the charge-density profile in "^®Ar is es¬ 
sentially related to the order of the proton states n\dj,i 2 and 
7r2ii/2, as well as the energy gap between the states. Among 
the selected Lagrangians, PKOl, PK03, PKDD and DD-ME2 
seem to support the emergence of the bubble-like structure 
in the proton-density profile of Ar with similar mechanism 
as shown in Fig. |3b). However, the DDRHF Lagriangian 
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TABLE I. Binding energies per particle EjA (in MeV) for N = 2& 
isotones from Z = 12 (Mg) to 20 (Ca) calculated by RH(F)B with 
PKAl, PKOl, PK03, PKDD and DD-ME2. The experimental data 
are extracted from Ref. ED and cr stands for the root-mean square 
deviations from the data. 



Exp 

PKAl 

PKOl 

PK03 

PKDD 

DD-ME2 

““Mg 

6.621 

6.625 

6.622 

6.524 

6.458 

6.456 

«Si 

7.416 

7.405 

7.392 

7.317 

7.301 

7.293 

44s 

7.996 

7.950 

7.966 

7.921 

7.896 

7.877 

46 Ar 

8.412 

8.366 

8.405 

8.379 

8.357 

8.332 


8.667 

8.674 

8.695 

8.683 

8.665 

8.644 

(T 


0.030 

0.021 

0.072 

0.103 

0.113 


PKAl, which has more complete meson-exchange diagram 
than the others, does not prefer the bubble-like structure to 
occur in the argon isotopes and such judgement is evidently 
supported by the fact that PKAl presents better agreement on 
the charge radii of argon isotopes with the data ll49l than the 
others [see Fig. [TJb)]. In particular, the proton shell Z = 16 
deduced from the systematical analysis of the experimental 
data (the beta-decay Q values, single proton separation en¬ 
ergies, and the energies of the hrst excited state) lISOll is prop¬ 
erly reproduced only by PKAl along the sulfur isotopic chain, 
from which the order of the states Ti\dii 2 and nlsi/x is decided 
certainly with the experimental evidence. That is, the proton 
state 7r2si/2 must be deeper bound than nld^ii to give the pro¬ 
ton shell Z = 16, and thus the occurrence of the bubble-like 
structure is prohibited. 



r(fm) r(fm) 


FIG. 5. (Color online) The charge-density distributions of N=28 iso¬ 
tones calculated by RHFB with PKAl, PKOl, PK03 and RHB with 
PKDD. 

To further probe the bubble-like structure, we extend the 
exploration from "^®Ar along the isotonic chain of A = 28, 
by taking the isotones from "^°Mg to ^^Ca. Table [n| shows the 



FIG. 6. (Color online) Canonical proton single-particle energies 
along the N = 2& isotonic chain calculated by RHFB with PKAl 
and RHB with PKDD. The lengths of thick bars correspond with the 
occupation probabilities of the proton orbits and the filled stars de¬ 
note the experimental data taken from Ref. ED. 


binding energy per particle E/A of the selected isotones calcu¬ 
lated by RH(F)B with PKAl, PKOl, PK03, PKDD and DD- 
ME2, in comparison with the data Il48ll . As seen from the root- 
mean square deviations in the last row, all the selected mod¬ 
els present appropriate agreement with the data, and PKAl 
and PKOl reproduce the binding energies better than the oth¬ 
ers. Turning to the charge-density profiles in Figs. |^a)- 
(d) which respectively show the results calculated by RHFB 
with PKAl, PKOl and PK03 and by RHB with PKDD, it 
is found that the calculations with PKOl, PK03 and PKDD 
present similar charge-density distributions with each other, in 
which the bubble-like structures, i.e., the central depressions, 
are predicted for the selected isotones only except '^^Ca [see 
Figs. I^b)-(d)]. While the central depressions exist only in the 
charge-density profiles oD^°Mg and “^^Si, as determined by the 
RHFB calculations with PKAl [see Fig. j^a)]. Such distinct 
deviations between the models can be understood from the 
proton single-particle spectra shown in Fig. It is seen that 
PKAl predicts deeper bound n2si/2 state than TT\dj ,/2 for the 
selected isotones and the state n 2 si /2 is gradually occupied by 
the valence protons since (Z = 16). As a result, the for¬ 
mation of the bubble-like structure is blocked in the isotones 
from to "^^Ca. Different from PKAl, the calculations with 
PKDD present an inversion on the order of the states n2s\i2 
and n\d 2 ,i 2 from Z = 12 ("^°Mg) to 14 ("^^Si) and then the va¬ 
lence protons are filled mainly in the orbit n\d 2 ,i 2 , leading to 
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the distinct central depressions in the charge-density profiles 
from to "^®Ar. Similar systematics are also found in the 
proton spectra determined by PKOl, PK03 and DD-ME2 as 
by PKDD. 


In fact, as shown in Figs. 2]|3 i.e., the charge-density pro¬ 
files and the corresponding proton spectra along the argon 
isotopic chain, the occurrence of the bubble-like structures is 
tightly related not only to the order of the states 7 r 2 si /2 and 
nldi /2 but also to the splitting between these two pseudo-spin 
partners. As shown in Fig. |^b), although the order of 7 r 2 si /2 
and ttI 6 / 3/2 is reversed at A = 22 , the emergence of the proton 
bubble-like structure in ""^“^"^Ar is still not favored very much 
[see Figs, [^b)-(d)] because of the fairly large occupations 
in 7 r 2 ii /2 induced by the pairing correlations, which is essen¬ 
tially influenced by the energy gap between the states. For 
the A = 28 isotones, if referring to the experimental data ED 
denoted by stars in Fig. the pseudo-spin symmetry related 
to the pseudo-spin partners 7 r 2 si /2 and nldi/i in "^^Ca is prop¬ 
erly reproduced by PKAl, whereas the calculations with the 
others show distinct discrepancy from the data, e.g., PKDD 
presents notable splitting between the states 7t2si/2 and nld 2/2 
as shown in Fig. [^b). 

In order to better understand the deviations between the 
models in predicting the occurrence of the bubble-like struc¬ 
ture, it is worthwhile to check the pseudo-spin symmetry in 
the relevant nuclei. From the Dirac equation, the single¬ 
particle energy of a state a can be expressed as. 


Fq. — Ek,a + Eq-q, -t Eqj^q -t Ep a + E^^q, -i- E^ a En at (13) 

where Fj. „ (0 = cr, u,p, tt. A) and F^ „ denote the contri¬ 
bution of the kinetic energy, potential energy and rearrange¬ 
ment terms, respectively. According to Fq. ( [T3| ) and using 
the canonical wave functions determined from the RH(F)B 
calculations, the contribution to the splittings of the pseudo¬ 
spin partners n 2 s\i 2 and nld^p are determined by the se¬ 
lected Lagrangians for '^^Ca. The results are shown in Ta¬ 
ble including the experimental values of the average en¬ 
ergy F = i2E„2st/2 + 4 F;jirf 3 ^J /6 and the splitting AF as a 
reference. Identical with the proton spectra shown in Fig. 

only PKAl properly reproduce the splittings between the 
pseudo-spin partners 7 r 2 si /2 and n\d 2 / 2 - Specifically, the term 
AEcr+oj extracted from the calculations with PKAl, namely 
the balance between the strong cr- and w-meson fields, plays 
an important role in reducing the pseudo-spin orbital split¬ 
ting, which indicates that PKAl presents a difference balance 
between the nuclear attraction and repulsion from the others 
ED. Similar systematics are also found in the detailed con¬ 
tributions of the splittings of the partners n 2 s\i 2 and nld 2/2 in 
'^^Ar. 

Combining the results of the argon isotopes and A = 28 iso¬ 
tones, it can be concluded that the bubble-like structure of the 
charge-density profiles is predicted to occur in the A = 28 iso¬ 
tones ^°Mg and “^^Si commonly by the selected RHF and RMF 
Lagrangians, due to the fact that in these two isotones both 
the proton states n2s\i2 and Ji\d 2 ,i 2 are not occupied. How¬ 
ever, for the popular candidate ^®Ar, the RHF model PKAl 
does not support the occurrence of the bubble-like structure in 
the charge-density profiles, and evidently it can provide better 


TABLE II. Average energy E = {2E,^2sii2 + 4Fii</3^2)/6 (in MeV) and 
contributions (in MeV) to the splittings AF of the pseudo-spin part¬ 
ner states 7t2si/2 and ;rld 3/2 in '**Ca, determined by the calculations 
with PKAl, PKOl, PK03, PKDD and DD-ME2. The experimental 
values of E and AF are extracted from Ref. ED 



Exp. 

PKAl 

PKOl 

PK03 

PKDD 

DD-ME2 

F 

-16.05 

-16.571 

-16.693 

-16.957 

-17.205 

-17.084 

AF 

0.36 

-0.170 

1.718 

1.716 

2.154 

1.516 


AF* 

- 

1.763 

1.152 

1.244 

1.229 

0.907 

AEp 

- 

1.366 

0.877 

0.640 

0.791 

0.489 

AF, 

- 

0.736 

0.569 

0.964 

- 

- 

AEa 

- 

0.149 

0.185 

0.197 

0.062 

-0.048 

AF« 

- 

-1.350 

-1.015 

-0.971 

-0.908 

-0.775 


- 

-2.961 

-0.185 

-0.492 

0.855 

0.828 


AF, 

- 

5.907 

9.656 

8.620 

24.838 

28.649 

AF^ 

- 

-8.868 

-9.841 

-9.113 

-23.983 

-27.821 


agreement with the data of binding energies and charge radii 
of the argon isotopes than the other RHF and RMF models 
(see Fig. [^, as well as the shell structures Z = 14 and 16 
nearby [see Fig. Qa)]. In addition, if starting from “^^Ca in 
which the pseudo-spin symmetry related to the partner states 
n 2 si /2 and tt 16 / 3/2 is demonstrated to be conserved experimen¬ 
tally, the neighbored Ar is expected to have nearly degener¬ 
ated pseudo-spin doublet 7t2si/2 and 7 : 16 / 3/2 and consequently 
the occurrence of the bubble-like structure will be blocked by 
the pairing effects which lead to the spreading of the valence 
protons over these two states. Among the selected models, 
only the RHF model PKAl presents consistent prediction on 
the conservation of pseudo-spin symmetry in ^^Ca and "^®Ar. 
On the other hand, the occurrence of the bubble-like struc¬ 
ture is also tightly related to the order of the states 7r2ii/2 and 
n\dj,j 2 - From Fig. Bb), some sulfur isotopes seem to have 
bubble-like structure, according to the proton configurations 
determined by the RMF model PKDD. While if referring to 
the existence of the Z = 16 proton shell as indicated by the 
experimental analysis m, the bubble-like structure will not 
allow to occur as well. Eventually, the nuclei "^°Mg and '^^Si 
are predicted to have the proton bubble-like structure not only 
from the existence of the distinct central depressions in the 
charge-densitity distributions, but also from the proton single¬ 
particle configurations. According to the half-life of '^^Si and 
"^‘’Mg respectively as 12.5ms and 170ns |[52l . "^^Si may be 
treated as a potential candidate of proton-bubble nucleus for 
experimentalists. 


IV. SUMMARY 

In this work we have studied the charge-density profiles and 
the proton spectra of the argon isotopes and A = 28 isotones 
with the relativistic Hartree-Fock-Bogliubov (RHFB) theory 
using the effective interactions PKAl, PKOl and PK03, and 
with the relativistic Hartree-Bogliubov (RHB) theory using 
PKDD and DD-ME2. It is found that both models can re¬ 
produce the binding energies and charge radii of the argon 
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isotopes with certain quantitative precision. Specifically, the 
PKAl effective interaction provides the best agreements with 
the data, particularly on the emergence of the proton shells 
Z = 14 and 16 nearby, and therefore the RHFB+PKAl model 
is supposed to be the most reliable one among the selected 
models. In the calculations with PKAl, the inversion of the 
proton orbits n 2 s\i 2 and n\d-ii 2 is not found in the argon 
isotopes to support the occurrence of the proton bubble-like 
structure. Along the isotonic chain of A = 28, fairly distinct 
central depressions are found in the charge-density profiles of 
^^Si and '^^Mg from the calculations of all the selected mod¬ 
els, which are mainly due to the fact that the proton orbits 
7r2ii/2 and n\d 2 ,j 2 is not occupied by the valence protons, and 
^^Si may be treated as a potential candidate of proton bub¬ 
ble nucleus with longer life time than "^°Mg. In addition it 


has been noted that another anti-bubble effect, namely the dy¬ 
namical correlation, would quench the bubble structure in the 
ground-state of ^'^Si strongly lfT9l l20ll . In present work, ^^Si 
is assumed to be the spherical nucleus, and perspectively it is 
interesting to test the existence of the bubble structure in "^^Si 
after taking the dynamical correlation into account. 
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